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Abstract. We present an analysis of eleven years of 
Stromgren by photometry of the red semiregular variable 
star AI Aurigae. An early period determination of 63.9 
days is confirmed by the long-term light curve behaviour. 
The light curve shows semi-regular changes with a mean 
period of 65 days reaching an amplitude of 0^6 in some 
cycles. The b — y colour changes perfectly parallel the V 
light curve, suggesting radial oscillation to be the main 
reason for the observed variations. We estimate the main 
characteristics of the star (mass, radius, effective tempera- 
ture) that suggest radial pulsation in fundamental or first 
overtone mode. 
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1. Introduction 

The light variation of semiregular variables (SR) of sub- 
type SRa and SRb is generally associated with the 
pulsation of these low and intermediate mass red gi- 
ants located on the asymptotic giant branch (AGB) of 
the Hertzsprung-Russell diagram. Since the characteris- 
tic timescale of the variations is between 20 and 2000 
days (usually hundreds of days), semiregulars are typi- 
cal targets of amateur observers estimating visual bright- 
nesses. There are very few high-precision light, colour, 
and/or radial- velocity measurement covering many cycles 
that confirm the pulsational origin of the variations. Very 
recently, Lebzelter and his co-workers (Lebzelter 1999ab, 
Lebzelter et al. 1999) presented high-precision photomet- 
ric and infrared radial- velocity data for a sample of bright 
SR stars. Their results both illustrate the ability of au- 
tomatic telescopes to monitor this type of variable stars 
and support the assumption of pulsation (see also Percy 
et al. 1996 on involving amateur photoelectric observers 
into regular observations of small-amplitude red variables, 
SARVs). 

However, one cannot ad hoc accept pulsation to be the 
main reason of variability. As has been noted by Percy 
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& Polano (1998), virtually every M giant is variable in 
brightness. This has been also supported by the large- 
scale photometric surveys (Hipparcos, MACHO, OGLE) 
leading to the discovery of thousands of M-type variables. 
It has been pointed out by, e.g. Lebzelter et al. (2000) 
and Kiss et al. (2000), that specifying the physical reason 
responsible for the variation requires simultaneous light, 
colour, and radial-velocity data, as even consecutive cy- 
cles may show significant deviations. This kind of long- 
term observation is crucial when studying the possible 
mechanisms affecting the light changes of SR variables. 
Besides pulsation, other contributors to variability cannot 
be readily excluded, such as time dependent surface inho- 
mogeneities due to large convective cells (Schwarzschild 
1975, Lebzelter et al. 2000), ellipsoidal deformation due 
to the presence of a close companion, or spots on a rela- 
tively rapidly-rotating star. Most recently, Koen & Laney 
(2000) gave a detailed list of considerations concerning the 
possible reasons of variability and noted, that there is still 
no definite conclusion on this issue. Until the application 
of these mechanisms is better understood, each case must 
be studied separately. 

AI Aurigae (HD 259342, 9"1 < V < 9^8, P = 63^9, 
spectral type M5III - Chuadze 1973) is a poorly-studied 
red semiregular (subtype SRa) variable star. In the Gen- 
eral Catalogue of Variable Stars (GCVS - Kholopov et al. 
1985-1988) a period of 63'?9 is listed with an early epoch of 
maximum (JD 2426029, due to Beyer (1937)). Since then 
there has been no other photometric measurement pub- 
lished for this star. Due to the relatively small range of 
light changes, it is also neglected by the international or- 
ganizations of amateur variable star observers (American 
Assocation of Variable Star Observers - AAVSO, Asso- 
ciation Francaise des Observateurs d'Etoiles Variables - 
AFOEV and Variable Star Observers' League in Japan 
- VSOLJ). The AAVSO International Database does not 
contain this star, AFOEV has 16 points which are wrongly 
identified with AI Aur (they are 2 mag brighter), and 
VSOLJ collected only 30 individual estimates between 
1986 and 1999. Similar neglect is present concerning other 
observations: there was only one radial velocity monitoring 
{{vr) = +63 km s~^. Feast et al. 1972), while Dickinson 
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& Dinger (1982) listed AI Aur among the negative de- 
tections during their H2O survey. There is no metallicity 
determination in the Uterature. Its infrared colours (see 
later) place the star among the "blue" SRs (Kerschbaum 
& Hron 1992) and the presently available data suggest the 
star to be a regular member of this group. 

The main aim of this paper is to present a continu- 
ous photometric monitoring of AI Aur between 1985 and 
1996 which reveals the associated colour variations in ad- 
dition to obtaining an accurate light curve. (We note, that 
the presented observations were acquired during a large- 
scale photometric survey of solar type stars (Lockwood 
et al. 1997), without any specific reason. Therefore, this 
study is a by-product of that long-term observing pro- 
gramme.) The observed behaviour of AI Aur is very likely 
to be mainly due to radial pulsation, most probably in fun- 
damental or low-degree overtone mode. The next section 
deals with the data aquisition, while the light- and colour- 
curves are discussed in Sect. 3 together with the physical 
parameters of the star. A summary of the conclusions is 
given in Sect. 4. 

2. Observations 

AI Aurigae was observed on 199 nights between 12 Oct 
1985 and 8 Dec 1996. The observations were carried out at 
Lowell Observatory with the 0.53 m telescope. The detec- 
tor was a single-channel photoelectric photometer (EMI 
6256S tube) equipped with Stromgren b and y filters. The 
measurements were obtained almost always through a 29" 
diaphragm; some measurc;mcnts in bad seeing were made 
through a 49" aperture, and in bright Moonlight a 19" 
aperture was occasionally used. 

The differential magnitudes were aquired relative to 
HD 46159 (l/=8™054, b - i/=(F623, G8III), while the 
photometric stability was checked against HD 259662 
(y=9™556, b - y=Cr639, KOIII). The data were corrected 
only for differential extinction using long-term monthly 
means for kb and ky (Lockwood & Thompson 1986). 
The standard transformations are represented with slight 
colour terms (roughly V = y + 0.05(6 — y)i, b — y = 
1.08(6 — 2/)j). The internal consistency is of order of 0™005 
as suggested by the rms errors of the means for the com- 
parison stars in y and 6 (ffi"0046 and 0^0048). The quoted 
standard values for the comparison stars were determined 
on four nights from standard star observations (extending 
to several very red Landolt stars). Due to the red colour 
of the variable, there might be some remaining systematic 
shift of 0""020-0""G25 in the V brightness. In any case, the 
scatter of the observations is essentially negligible com- 
pared to the range of the variations. The data (mean val- 
ues of var— comp and var— check) are presented in Table 
1 (available electronically and/or in printed version?) and 
plotted in Fig. 1. 
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Fig. 1. The light and colour curves of AI Aur between 
1985 and 1996. Note the 600-day break between MJD 
48000 and 48600 



3. Discussion 

In the following, we characterize the observed light varia- 
tions, discuss the importance of the colour measurements 
and compare AI Aur with other semiregulars with avail- 
able similar data series. 



3.1. Period analysis 

Data presented in Table 1 are plotted in Fig. 1. Note that 
some parts are excluded from the presentation due to the 
sparse light curve coverage (between MJD 48000-48600). 
However, when analysing the whole dataset, we have also 
used those data that are not presented in Fig. 1. Qualita- 
tively the light curve is that of a typical semiregular with- 
out strictly repeating cycles. This has been investigated in 
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detail by means of standard Fourier analysis implemented 
with Period98 of Sperl (1998). 

First, we have calculated the power spectrum of the 
whole dataset. It is shown in Fig. 2, where the insert 
shows the window function. One main peak is present 



at f=G. 01604 c/d (A=0™14) corresponding to a period of 
62'?345. After a simple prewhitening with this frequency, 
the power spectrum of the residuals contains a closely sep- 
arated frequency (f=0. 01548 c/d, A=0^10). However, it 
is a direct consequence of the unstable period, since the 
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Fig. 2. The Fourier spectrum of the whole datasct. The 
insert shows the window function with the same frequency 
scale 



whole dataset could be fitted neither with a single period 
nor with a sum of two close periods. Therefore we stud- 
ied the light curve in four separate subsets (each being 
about 1000 days long) plotted in the left panels of Fig. 3. 
We performed similar frequency analysis as for the whole 
light curve and the results can be summarized as follows: 

1. There is a well-defined main peak in every Fourier 
spectra representing the average cycle length varying 
around the mean period of ~65 days. The fitted curves 
(solid lines in Fig. 3) have the following periods and 
amplitudes: 1 - 63'?2, a"2; 2 - 67^1, ffi"13; 3 - 100 d, 
0^12 (the few data points may result in wrong result); 
4 - 65 d, 0™15. The typical period uncertainty is about 
1 d (except for the third subset). 

2. The first subset can be fitted very well with two har- 
monic components, however the second one is simply 
due to the variation of the mean brightness. Later data 
do not support the presence of this false period (~860 
days) . 

3. The phase diagrams of the individual subsets (right 
panels of Fig. 3) suggest the average light curve shape 
to be quite stable. The data were phased with those 
periods listed above. 

4. The prewhitened spectra of the subsets do not contain 
any other significant peak higher than 0™08. We have 
not tried to use the low- amplitude component (s), be- 
cause the residual light curves show irregular rather 
than cyclic behaviour. We feel that when analysing 
light curves of SR variables the noise level in the fre- 
quency spectrum is only partly due to the observa- 
tional errors and it would be misleading to accept ev- 
ery peaks in the spectrum being much higher than 
the mean error of the observations (about O^'Ol in our 
case). The star itself may grossly increase the noise of 
the light curve through its erratic behaviour caused, 



e.g., by the mild chaos affecting the pulsation (see e.g. 
Buchler & Kollath 2000). It is well established, e.g. 
for Mira stars, that radial pulsation alone may pro- 
duce irregularities in the light and/or radial velocity 
curves (for instance, through shock wave propagation 
in the extended atmosphere, Bertschinger & Cheva- 
her 1985, Fox & Wood 1985). The complicated stel- 
lar responses on other physical mechanisms (coupling 
between pulsation and convection, possible departures 
from the spherical symmetry) prevent drawing a sim- 
ple and straightforward conclusion on the irregulari- 
ties. Most recently, Lebzelter et al. (2000) revived the 
old idea proposed first by Schwarzschild (1975) on the 
large convective cells introducing some irregularity in 
the light variation. However, in the case of AI Aur, 
only speculations can be drawn. 
5. The results slightly depend on the selection of the sub- 
sets (we have tried different lengths starting from one 
year to the finally accepted 1000 days); on the other 
hand, they are strongly affected by the lack of points 
in the second half of the data. 

Considering these points we conclude that the light 
curve cannot be described with simple sums of harmonic 
functions, but it has a slowly and apparently irregularly 
variable characteristic cycle length (instead of period) near 
65 days. The period analysis of the 6-filter data leads to 
the same conclusion. 

3.2. Colour variations 

The conclusions presented so far are typical for semireg- 
ular variables (see, e.g. Kiss et al. 1999). The novelty in 
our results is the presence of colour measurements. To our 
knowledge, there have not been really long-term (covering 
tens of cycles) time-series multicolour data in the litera- 
ture for any semiregular variable star. Although Hipparcos 
Epoch Photometry (ESA 1997) contains typically 100-200 
points for thousands of variables, the limited time span 
does not permit accurate period determination for red 
variables. High-precision V-band observations of different 
SRs were presented, for instance, by Percy et al. 1989, 
Cristian et al. (1995), Percy et al. (1996) and Lebzelter 
(1999), but none of these studies dealt with the colour 
variations. 

One of the most extensive time-resolved multicolour 
photometric surveys was published by Smak (1964), who 
reported U BV observations of twenty-nine Mira stars and 
twelve scmircgulars during two years, covering typically 
only one cycle per star. He did not find any clear corre- 
lation of the B — V colour and the V brightness for the 
SR stars. Jerzykiewicz (1984) studied the light and colour 
variations of HD 157010 (V818 Her), concluding that the 
star's UBV magnitudes are effective temperature param- 
eters through the sensitivity of the continuous spectrum 
and TiO blanketing. Wisse (1981) presented UBV obser- 
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Fig. 3. The four subsets and the corresponding phase diagrams of AI Aur 



vations for 35 SRs for classification purposes. He found 

that well-defined correlations and anti-correlations exist 
among the variables, where the phase difference of the V 
and B — y maxima is either or 0.5. Cadmus et al. (1991), 
in their study of three SRs, noted only that the colour 
changes were small compared to the changes in brightness. 
Narrow-band observations of SR variables in the literature 
are quite scanty. An interesting application can be found 
in Wasatonic & Guinan (1998), who used three colors of 
the Wing near-infrared system to trace the temperature 
and radius variations of V CVn. These authors measured 
the highest temperature (from an infrared index) around 
the maximum light. 

The most important feature of our observations is the 
almost perfectly parallel light and colour curves. This is 
a characteristic behaviour in radially pulsating variables, 
where the temperature maximum occurs close to the max- 
imum light. It means that the star appears bluer when it is 
at maximum light, which is in good agreement what is ex- 
pected from pure thermal variations driven by pulsation. 
In the case of FS Comae (SR star with PRi58 d), Tor- 
res et al. (1993) reached the opposite conclusion based on 
U BV data series. The star showed reversed colour changes 
explained by the fact that B — V colour in FS Com is a 
molecular absorption index, rather than a temperature in- 
dicator. The Stromgren h filter (centered near 4700A) has 
a relatively clear spectral region in the center of the pass- 



band, but the wings on both sides have moderately-weak 
TiO absorption, particularly the 2.0 transition with band- 
head at 4761A. The y filter is not so strongly affected, with 
only the usual Fe, Ti, etc. lines. Unlike the B — V colour, 
& — y is completely unaffected by the wide TiO bandhead 
at ~5167A. We conclude that the observations imply ra- 
dial oscillation with dominant thermal effects as the main 
reason for the light variability. The ^^0™! amplitude of the 
h — y curve corresponds to a temperature change of <300 
K (as estimated from the synthetic colour grids of Kurucz 
1993). 



3.3. Stellar parameters 

Adopting radial pulsation for AI Aur, further constraints 

can be drawn on the basic stellar parameters. Unfortu- 
nately, there is no parallax measurement for AI Aur, thus 
only spectroscopic parallax or various empirical period- 
colour-luminosity relations (e.g. Barthcs et al. 1999) can 
be used to determine its luminosity. We have chosen the 
former approach, because we wanted to estimate stellar 
properties without any assumption on its pulsation. In 
the following we have neglected the interstellar reddening 
since the observed mean b — y colour (1™35) is close to 
the expected value for an M5 giant star (see Jorissen et 
al. 1995 for Stromgren photometry of red giants). 
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The IRAS [12]- [25] colour is -1""41. Following the def- 
initions by Kerschbaum & Hron (1992), this means that 
AI Aur belongs to the "blue" semiregulars (no indica- 
tion for circumstcUar shells, P<150'^, T^s > 3200 K). 
The infrared brightnesses {K=2'^71, /=6™18), taken from 
the IRC catalogue (Neugebauer & Leighton 1969), com- 
bined with the mean V magnitude imply colour indices 

V - K=&^8 and V - J=3^3. The latter value results 
in a reff=3420 K using the temperature scale by Dumm 
& Schild (1998). A less biased estimate can be inferred 
from twenty M5III-type stars in the sample of Dumm & 
Schild (1998): (Teff(M5777))=3520±100 K. For an inde- 
pendent check we have used the {V — K)-{710 — 888) and 
(/ - ii:)-(710 - 888) relations in Alvarez et al. (2000) re- 
sulting in refF=3300 K from their TefF-(710 - 888) cali- 
bration. Similarly, one can get a mean radius and mass 
for the given spectral type of (R) = 123 ± 34i?0 and 
(M) = 2.0 ± O.8M0. Adopting these spectral type-mass 
and spectral typc-radius values we calculated a pulsational 
constant Q{= P ■ ^p/pq)=0.0651^;^^. Recently, Percy 
& Parkes (1998) discussed the pulsation modes in small- 
amplitude red variable stars reaching a conclusion that 
some stars are likely to pulsate in up to the third overtone 
mode. Within these frameworks AI Aur seems to pulsate 
in the fundamental or first overtone mode as is suggested 
by theoretical models of Xiong et al. (1998) or Ostlie & 
Cox (1986). The overtone pulsation is favoured by the po- 
sition of AI Aur in the ii'-band P-L diagram of long-period 
variables derived by Bedding & Zijlstra (1998) from Hip- 
parcos parallaxes. For this, we estimated Mk from the 

V — K colour and My-spectral-type calibration of The et 
al. (1990). The resulting Mr = -7^5 and log Pwl.8 place 
AI Aur close to the upper sequence in Fig. 1 of Bedding 
& Zijlstra (1998), which may be interpreted as a conse- 
quence of pulsation in a different mode than that valid 
for most Mira stars (represented by the lower sequence in 
Fig. 1 of Bedding & Zijlstra 1998). 

Finally, there is an interesting period-gravity relation 
for a wide range of radially pulsating variable stars pre- 
sented by Fernie (1995), which can also be used to test 
the assumption of radial oscillation. The adopted mass 
and radius give a log g(=0.56±0.4 which is coincidentally 
the same as predicted by Eq. 1 of Fernie (1995). Plotting 
AI Aur in Fig. 1 of Fernie (1995), its position is as de- 
viant as that of Mira itself suggesting the first overtone to 
be somewhat more likely. However, the universality of this 
period-gravity relation has no firm theoretical background 
and verification; therefore, this comparison should be con- 
sidered only as a possible hint for the mode of pulsation. 

The presented considerations are on the whole consis- 
tent with the recent observational and theoretical results 
regarding the mode of pulsation in Mira and semiregular 
variables. For example, Feast (1996) found that semireg- 
ular variables, independently of their metallicity, pulsate 
probably in the first overtone mode. Further supporting 
arguments were listed by Feast (1999). However, semiregu- 



lars form a quite heterogeneous group, in which stars may 
pulsate in the fundamental, or 1st, 2nd or even 3rd over- 
tone mode, as has been clearly demonstrated by Wood et 
al. (1999). That is why every individual case study has to 
be performed without any definite preconception. 

4. Conclusions 

The main results presented in this paper can be summa- 
rized as follows: 

1. New Stromgrcn hy photometric observations of the 
semiregular variable AI Aurigae are presented and dis- 
cussed. The period analysis of the whole light curve con- 
firms the early period determination of 63.9 days. How- 
ever, the light curve shows typical semiregular behaviour 
with changing amplitude and cycle length. 

2. The b — y colour measurements revealed highly paralel 
light and colour variations. Since the h — y colour is consid- 
erably good temperature indicator (rather than molecule 
indicator, as is the B — V colour), the phase dependent 
light and temperature changes resemble classical radially 
pulsating stars. The observations imply radial oscillation 
with dominant thermal effects as the main reason for light 
variability. 

3. We estimated the basic stellar parameters, such as the 
effective temperature (3520±100 K), mean radius (123±34 
R0), mass (2.0±0.8 Mq). The spectral type-mass amd 
spectral type-radius values give a pulsational constant of 
0.065j;Q Q3. The results are consistent with radial pulsation 
in the fundamental or (more likely) first overtone mode. 
The period-gravity relation of Fernie (1995) also supports 
this conclusion. 
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